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ABSTRACT 

We discuss detailed fits of the BATSE and PVO gamma-ray burst peak-flux distributions 
with Friedman models taking into account possible density evolution and standard candle 
or power law luminosity functions. A chi-square analysis is used to estimate the goodness of 
the fits and we derive the significance level of limits on the density evolution and luminosity 
function parameters. Cosmological models provide a good fit over a range of parameter 
space which is physically reasonable. 
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1. INTRODUCTION 

Gamma ray burst sources are distributed with a very high level of isotropy (Fishman et al. 

, 1994), which is compatible with either a cosmological origin or an extended galactic halo 
origin. The brightness distribution is another indicator used to characterize the spatial 
distribution in distance, and this can be used to further test the distance scale hypotheses. 
This is generally done by investigating the functional behavior of the integral number N of 
sources with peak photon flux rates P above a certain value, N{> P), or of the peak count 
rate divided by the threshold rate N{> Cmax/Cmin)^ or the corresponding differential 
distributions. Comparisons of observed versus expected values in Friedman cosmologies 
have been discussed, e.g., by Mao & Paczyhski (1992), Dermer (1992), Piran (1992) 
and Wasserman (1992). Statistical fits to a log iV — logP or log iV — log C distribution 
have been done by Loredo & Wasserman (1992), Wickramasinghe, et al. . (1993), Cohen 
& Piran (1994), Emslie & Horack (1994), Horack, Emslie & Hartmann (1995), Fenimore 
& Bloom (1995), Rutledge, et al. (1995). One of the main questions that such fits must 
address is the size of the parameter space region which is compatible with a cosmological 
distribution, and whether such parameters are reasonable. If the acceptable region contains 
physically plausible parameters and is not too restricted, one may assume the consistency 
of the observations with a general type of models; if on the other hand the acceptable 
region is very small and/or populated mainly by physically implausible parameters, fine- 
tuning would be required to fit the observations, and the case for consistency with those 
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models is weaker. Such consistency, and absence of fine-tuning, is a requirement expected 
of any successful model of the GRB distribution, whether cosmological or galactic. Here, 
we shall address only the question of the consistency of the number distribution under the 
hypothesis of a cosmological distribution. 

So far, most cosmological fits have been made with relatively specialized models, 
generally either with non-evolving or evolving density standard candle models, or with 
non-evolving luminosity functions. Limits on the luminosity function were investigated in 
cosmology with a pure density evolution by Horack, Emslie & Hartmann (1995) using a 
method of moments. In Euclidean space, limits have been investigated by Horack, Emslie 
& Meegan (1995), Ulmers & Wijers (1995) and Ulmer, Wijers & Eenimore (1995). Most 
cosmological calculations have used either the IB or the 2B BATSE data base, and did 
not include the PVO information (see, however, Cohen & Piran, 1995, Eenimore & Bloom, 
1995). In the present paper we make detailed chi-squared fits of the observed brightness 
distribution directly to specific models of the cosmological burst brightness distribution. 
We use both standard candle and power-law luminosity function models with a density 
evolving as a power law of the scale factor, for a wide range of density evolution exponents, 
luminosities and luminosity spreads, assuming either a brightness limited or redshift limited 
cases for various maximum redshifts for the source distribution. This is done both using 
the BATSE 2B catalogue of sources (Meegan, et al. , 1994), and combining the BATSE 
catalogue with information published for the PVO counts (Eenimore and Bloom, 1995). 
The significance levels of the various cosmological fits is discussed for both the 2B and the 
expanded burst sample. 


2. COSMOLOGICAL DISTRIBUTION MODELS 


Analytical expressions for the integral burst number counts N{> P) with peak 
photon flux rate in excess of P (units of photon cm“^ s~^) were discussed by Meszaros & 
Meszaros , 1995 (MM95) for arbitrary Eriedman models with zero cosmological constant 
(in that paper C was used for the photon flux, but here instead we use P for the photon 
flux to avoid confusion with the more common usage of C [ s“^] as the count rate). As 
discussed in MM95, effects of a non-flat cosmology (Oq < 1) are small, and to a first 
approximation can be neglected. Below we assume Oq = 1 everywhere. The effect of a 
pure density evolution is approximated through a dependence 

n{z) = no{l + z)^ , (1) 

where n is the physical burst density rate in cm“^ Uq is the density rate at 2 ; = 0 

(D = 3 corresponds therefore to a non-evolving, constant comoving density). Eor a source 

emitting C photons s“^ with a power law photon number spectrum Ci^ oc (i.e. 

a = 0 corresponds to a flat power-per-decade spectrum), assuming most of the photons 
are collected in an energy range where a ~ constant a K-correction is necessary (e.g. Mao 
& Paczyhski , 1993, Dermer, 1992). This can be folded in with the density evolution by 
using an effective scale factor exponent D^ff = {D -|-a — 1) (MM95) and this K-correction 
is small or does not apply to most bursts, for which a ~ 1 in the range 50-300 keV where 
BATSE collects most of the GRB photons used to determine the peak flux P (e.g. Band, 
et al. , 1993). The photon luminosity function in the 50-300 keV range is represented by 
either one of the two forms. 


$(£) = 


noS{C — Co) , (standard candle) ; 

nCl^^^iCICrmn)~^ for Cmin < C < Cmax (pOWer law) 


( 2 ) 
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For the rest of the paper we take C to be the burst peak photon luminosity [ s“^], P to 
be peak photon flux [ cm“^ Uq is the physical density of bursts per year at 2 ; = 0, 

fi = no(l — l3)/(K^~^ — 1) and K = Cmax/^min gives the effective spread of the intrinsic 
luminosity function in the power law case. 


The integral number distribution of bursts per year with peak flux rate above P 
is given by (MM95) as 


N{> P) 


AtT 

3 (47rP)3/2 


(3) 


where Ce = Co^ = Ug for the standard candle, Cg = Cmin^ = n for the power 

law luminosity function, and / is a dimensionless analytical function of S and the model 
parameters, i.e. the luminosity function parameters Cg or K^Cmim the density 

evolution parameter D. In the redshift limited case, the maximum source redshift z^ax is 
an additional parameter, and the expression corresponding to equation (3) is given in the 
appendix of MM95. The differential number distribution N(P) = —dN[> P)ldP can be 
obtained from the integral expressions through differentiation. 


3. STATISTICAL MODEL FITS 

For the numerical fits we used the BATSE 2B catalog available electronically. 
The 1024 ms peak fluxes P (photons cm“^ s“^) were used, and only events with peak 
count rates divided by threshold count rates CmaxICmin > 1 were included, where Cmin 
is the published count threshold for each event. The 2B sample with this criterion consists 
of 278 entries in the catalog. Applying the efficiency tables published with the catalog to 
correct for detector inefficiency near the trigger threshold, the nominal number of bursts 
accumulated by BATSE over a period of two years with peak fluxes above logP > —0.6 is 
369. We chose for these bursts a binning equidistant in log^^Q P, with step size 0.2 between 
-0.6 and 1.2, which gives 9 equal bins with a minimum number of 7 events per bin (in the 
highest S bin, log^^Q P = 1.0 to 1.2) for the two year 2B sample. The fits were made to 
the differential burst number distribution N(P) as a function of peak photon flux P (since 
only in the differential distribution may the bins be considered independent of each other 
for a fit) and the errors in each bin were taken to be the square root of the number of 
events in that bin. 

Some of the fits were made using an extended 2B plus PVO sample. Eor the 
PVO events, we used the PVO portion of Table 2 of Eenimore and Bloom (1995) [EB95] 
for log;^Q P > 1.2. A number of subtle issues concerning a matching between the different 
PVO and BATSE data sets are discussed by Eenimore, et al. , 1993, who indicate that 
systematic uncertainties of ±10% in the relative normalization cannot be ruled out. The 
matching of the level of the BATSE and PVO curves was taken directly from EB95. The 
PVO data was rebinned, ignoring PVO bursts below log^^Q P = 1.2 so as not to count 
twice, and its level was renormalized so that the matching 2B data had the same level as 
in the original 2B catalog, i.e about 2 years. The errors for the PVO sample were also 
renormalized taking into account the fact that data had accumulated over more than ten 
years in the PVO case, keeping the relative errors the same. We used 5 bins in the PVO 
range, so that the combined 2B±PVO fits have 9±5=14 bins, reaching up to log^^Q P = 3.0. 

i) The SC fits (standard candle with density evolution) involve the fewest pa¬ 
rameters: the photon luminosity Cg (ph/s), the density Ug and the density evolution 
parameter P, under the brightness-limited assumption. Eor the 2B sample between peak 
fluxes —0.6 < log P < 1.2 the free parameters are p = 3, the degrees of freedom are 
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/ = 6 and the best Xred (reduced chi-square or divided by degrees of freedom) is 0.85 
at the innermost mark. The Icr, 2(7, 3(7 significance contours were determined using the 
standard prescription (e.g., Press, et al. , 1986, or Lampton, Bowyer and Margon, 1976). 
The fit (Figure la) is good over an elongated region describing a relation between the 
luminosity and the density evolution. For faster density evolution n oc (1 -|- z)^ (larger 
D) the luminosity must increase because most sources are at larger redshifts, while for 
slower or negative evolution the luminosity must decrease, since most sources are at small 
redshift {D = 3 is constant comoving density). The optimal fit is obtained for D = 3.5 and 
Co ~ 10^^ s“^. This luminosity is close to the SC value deduced, e.g. by MM95 and some 
previous authors as well (corresponding to Lq ~ 10^^ erg s“^ for typical photon energies of 
0.5 MeV). However the preference for D = 3.5 was not, as far as we can tell, encountered 
in previous fits. The Xred ^-round the best fit minimum is 0.85; however, the Icr region 
around it is rather large, even if not very wide, so this preference is not strong. 

For the SC fits using the 2B-|-PVO sample, p = 3, / = 11 and the fits are shown 
in Fig. lb, with a best Xred ~ 0.62 at the central mark enclosed by its Icr, 2a, 3a contours. 
We note that this ignores any possible systematic errors in matching BATSE and PVO 
beyond what is done in Fenimore et al. , 1993, and Fenimore and Bloom, 1995. If any 
extra errors were present, they could in principle increase the size of the significance regions 
discussed below (e.g. it might add an extra free parameter for the relative normalization). 
However, such errors are extremely difficult to quantify without going into additional 
details of the instruments, and we follow Fenimore and Bloom (1995) in adopting their 
relative normalization as adequate without further manipulation. The effect of the rare 
high flux PVO bursts satisfying a tight N oc correlation at 1.2 ^ log;^o P ^ 3.0 

is to improve the best fit (lower X^red) place it at a somewhat smaller luminosity 

£0 ~ 5 X 10^® s“^ and closer to comoving constant density evolution, 17 ~ 3. This is 
in good agreement with Fenimore and Bloom’s (1995) value of To ~ 5 X 10^** erg s“^. 
However, as seen from Fig. lb, the Icr region around this best fit minimum is compatible 
with both larger and smaller Co and D. In contrast to the pure 2B fit, however, the joint 
1(7 upper limit for Co and D are Co ^ 5 X 10^^ s“^, D 4.5 (or 3a joint upper limits 
£0 ^ 5 X 10^® s-i, D ^ 5). 

ii) The PL fits (Power-Law luminosity function bounded between Cmin and Cmax 
and including density evolution) are shown in Figures 2 and 3 in the brightness limited 
case. The parameters of the fits are Cmax^ 77, Uq, D, where K = Cmax jCmin-, and for the 
2B sample p = 4, / = 5. If the index /3 is taken as an additional variable parameter, the 
fits maximize at slopes significantly steeper than -5/2, and this results in Cmin dominating 
the luminosity function over the whole range of S, giving essentially a standard candle 
case. As discussed in MM95 (see also Ulmer and Wijers, 1995, Wasserman, 1992), this is 
because a luminosity function slope /3 = —5/2 reproduces directly the Euclidean integral 
distribution slope -3/2 (which is a differential distribution of slope -5/2). However, the 
quality of the fits with such very steep slopes (or SC cases) is not significantly different 
from those for a fixed slope of /3 ~ 1.88, in the sense that in the brightness limited case 
both give Xred ^ index close to -1.88 is suggested for a power law luminosity slope, 

since the slope of the integral distribution at low S is approximately -0.88 (e.g. Meegan 
et al. , 1992; see also Wasserman, 1992). An illustration of this is given, e.g. in Eig. 1 of 
MM95. Either steeper or shallower slopes would lead to an effective standard candle case 
below Pturn ~ 6 cm“^ s“^. Since there is no way to either rule out or prefer the physically 
interesting slope of -1.88 in the brightness limited case, we discuss here fits which take a 
fixed value /3 = 15/8 ~ 1.88. Given such an index the interesting question is what can be 
said about the lower and upper limits Cmin or Cmax^ or Cmax and the intrinsic spread K, 
and how do the fits compare with SC fits or K = 1 fits of the same slope. The 2B PL 
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fits (Figure 2) in the Cmax^K parameter plane show a broad inverse correlation between 
the allowed values of Umax and K. For increasing values of the density evolution index 
D the required values of Cmax also increase, as one would infer from the previous SC fits, 
since the sources are more distant and need to be more luminous. The best Xred 
cases D = 2,3,4 (in panels 2a,2b,2c) are ~ 0.82 in all three cases, at the (varying) location 
of the innermost mark surrounded by the Icr, 2(7, 3(7 contours. Fits have also be obtained 
for other values of —2 < D < 5, the trend being apparent from the three values shown 
here. These 2B best fits are obtained for values of Cmax^K ~ 2 X 10^^, 30 for D = 2, 
Cmax 7 K ~ 3 X 10^^, 15 for = 3, Cmax-,K ~ 8 X 10^^, 5 for D = 3. However, the 2B 
fits are not very strongly constrained in the Euclidean -3/2 region, and the joint Icr limits 
could extend to relatively large values of Cmax and of the intrinsic luminosity ratio K. 

The similar PL fits for the 2B-|-PVO sample in the brightness-limited case, with 
p = A, f = 10, are shown in Figs. 3a,b,c for D = 2,3,4. The presence of the bright 
PVO bursts constrains the high luminosity portion of the fits, and brings the best fit 
minima towards somewhat lower luminosities Cmax and lower luminosity ratios K. The 
best fit joint minima are Xred ~ 0.62,0.62,0.70 near Cmax,K ~ 7 X 10^®, 10 for D = 2, 
CmaxiK ~ 8 X 10^®, 3 for D = 3, Cmax-,K ~ 1.5 X 10^^, 1 for D = A. While the best fit 
values of K are relatively small, it is to be noted that for all three D the joint Icr upper 
limits allow intrinsic luminosity ratios of order K ~ 10^, while the 3(7 upper limits are of 
order K lOh 

hi) The PLZ fits (luminosity function with density evolution and redshift cutoff 
model) differ from the above in that a maximum redshift Zmax is included as a parameter 
in the expression for the distribution N (in the integral I of equation [3], see MM95). The 
parameters are Cmax^ 77, Uq, 77, Zmax^ so p = 5, and for the BATSE 2B sample (/ = 4) the 
fits are shown in Eigs. 4a through 4f, for D=2,3,4 and two particular choices of Zmax = 2, 6. 
The best Xred values marked inside the Icr, 2(7, 3(7 contours are 1.1, 1.0, 1.2 for 17 = 2,3,4 
{zmax = 2) and 1.2, 1.2, 0.8 for 17 = 2,3,4 {zmax = 6). In the Cmax^K plane the results 
are shown in Eigure 4a,b,c (left panels) for the case Zmax = 2.0, while Eigure 4d,e,f (right 
panels) is for Zmax = 6.0. Three values of the evolution index D are shown. In the D = 2 
case most of the sources are nearby, so the 2B fits with Zmax = 2 (Eigure 4a) are not 
very different from those for Zmax = oo (Eigure 2a) except for somewhat lower Cmax^K 
values. Eor 17 = 4 however, most sources would be farther and the Zmax = 2 case (Eigure 
4c) restricts the low luminosities and requires higher Cmax^K values than the brightness- 
limited Zmax = oo case of Eigure 2c. Eor Zmax = 6, however, the results (Eigure 4 d,e,f) 
are fairly close to the brightness limited case (2a,b,c), except for 17 = 4 where there is 
some restriction at low luminosities. 

Eor the 2B-|-PVO sample, the PLZ power-law maximum redshift fits provide some 
additional restrictions in the Zmax = 2 case. The PLZ fits are shown in Eigs. 5a through 
f, where p = 5, f = 9, and the Xred values are 0.7, 0.7, 0.9 for Zmax = 2, and 0.7, 0.8, 0.7 
for Zmax = 6. While 17 = 2 is fairly similar for both redshifts (sources are close by in both 
due to evolution) the 17 = 4 case (sources far due to evolution) is constrained by the PVO 
data to have larger K and Cmax values, and gives a joint 3(7 lower limit of K ^ 5 for the 
Zmax = 2 case. For the Zmax = 6 case, however, the 2B-|-PVO fits (Figure 5d,e,f) are not 
very different from those in the brightness limited case (Figure 3), except for a not too 
significant preference towards somewhat larger Cmax^K- 
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4. DISCUSSION 


The fits presented above show that a cosmological interpretation is compatible 
with the data under a variety of assumptions. Good fits to the observed differential dis¬ 
tribution of bursts N(P) as a function of peak photon flux P are obtained both under 
a standard candle (SC) and under a power-law (PL) luminosity function assumption, as 
well as for redshift-limited power-law luminosity function (PLZ). Fits were obtained for 
a range of density evolution indices D, defined through a physical density dependence 
Uo Gc (1 z)^ where D = 3 is equivalent to a non-evolving, constant comoving density 

case (i.e. S = D — 3 = 0, where S is the index for the comoving density ricom = ?^o(l + 
sometimes used in the literature). The BATSE 2B fits, which have many weak bursts but 
few very bright bursts, do not constrain D except through an inverse functional dependence 
on the luminosity. For the 2B-|-PV0 fits, which include a number of PVO bright bursts, a 
3a upper limit to the SC luminosity and the evolution index are obtained, jCq ^ 5x 10^® s“^, 
D b, "With optimal values at P '0 b 'X. fO s , P) 3. 

In the 2B and 2B-|-PVO samples there is no conclusive evidence for a lower 
cutoff Pmin (or a flattening of the integral distribution), the uncertainty being due to 
large and uncertain trigger corrections near threshold. In the brightness-limited case, the 
lowest peak flux P used also corresponds to the largest redshift observed, via the relation 
Co = 47ri?g[(l + — VyP = 4 X 10^’^/i“^[(l -|- x)^/^ — VyP s“^, where Ro = 2c/iLo is 

the Hubble radius, and Ho = 100 h km/ s/ Mpc is the present Hubble constant. Thus by 
re-plotting Figure 1 as a function of z{Pmin) and taking into account differences in notation 
for D one can also obtain a plot similar to the figure 3 of Cohen and Piran (1995), although 
the sample criteria and cuts are somewhat different. Our best fit Co varies with D (Figures 
la,b), and log^^Q = —0.6 corresponds in Figure lb to z^ax ^ 1-9 (10.; 0.8) for our 
best fit (ilci) values of Co ^ 5 X 10^® s“^, D ~ 3 (Co ~ 5 x 10^^ s“^, D ~ 4.5; Co ~ 
1 X 10^® s“^, D ~ 0.0). (To be specific, here we arbitrarily took a lower value of D ~ 0 
but from Fig. lb one sees that the lower limit could actually be smaller and could extend 
to D ^ —2 and values of Co ^ 10^® s~^)- Thus, if an intrinsic energy-stretching of the 
time profiles exists indicating a maximum SC redshift z^ax ~ 6 (as argued by Fenimore 
and Bloom, 1995), this could be easily accommodated within our 2B-|-PVO Icr SC limits 
with a density evolution faster than comoving constant, D >3. 

The fits with a power law (PL) luminosity function (equation [2]) bounded be¬ 
tween Cmin and Cmax behave in a manner which is qualitatively similar to that in the SC 
case. The fits presented are for a luminosity function index /3 = 15/8. In the brightness 
limited case using the BATSE 2B sample, there is an inverse correlation between max¬ 
imum luminosity and the evolution index D. (As discussed in Meszaros and Meszaros 
, 1995 [MM95], for the Euclidean P~^l'^ part of the integral distribution the behavior is 
dominated by the large luminosity sources C ~ Cmax^ if fi < 5/2). While for 2B there is 
a preference for a ratio of intrinsic luminosities K = Cmax jCmin of order 10-30, the Icr 
upper limits are compatible with much higher values. However, using the PVO data as 
well, we obtain more specific constraints on K and Cmax- For D = 2, 3, 4 the best fit 
K and its Icr upper limits are (5,100); (2,60); (1,30) (see Figure 3). These results for 
D = 3 (nonevolving density) are in significant agreement with the results obtained from the 
method of moments on the observed luminosity distribution by Emslie & Horack (1994), 
Horack, Emslie & Meegan (1995) and Horack, Emslie and Hartmann (1995), using the 2B 
sample. They are also compatible with the Euclidean distribution results of Hakkila, et al. 
(1994) using 2B, as well as results on the observed and/or intrinsic luminosity distributions 
by Ulmer & Wijers (1995) using 2B, Ulmer, Wijers & Fenimore (1995) using 2B-|-PVO 
data and Hakkila, eA al. (1995) using 3B -|- PVO data. It is worth noting that while some 
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of these groups have put limits of a factor ~ 10 for the width containing 90% of the bursts 
in the observed luminosity distribution, our widths here refer to the intrinsic luminosity 
distribution. Also, note that our width limits apply for the particular luminosity function 
slope /3 ~ 15/8 which fits the low flux end of the observed number distribution. As empha¬ 
sized in MM95, for such a slope each luminosity decade contains only about 10% as many 
bursts as the previous lower luminosity decade, so ~ 90% of the bursts are automatically 
in the lowest decade. For other slopes /3<lor>5/2 the intrinsic widths could be much 
larger, since the upper or lower ends of the luminosity function dominate and one is dealing 
effectively with standard candles. 

By assuming a redshift cutoff to the source distribution, the corresponding power- 
law luminosity function fits (PLZ) involve more than just a simple relation relation between 
Pmin and Zfnax- As discussed in MM95, even for P > Pmin the integrations over the 
luminosity functions depend on z^ax in a nontrivial manner. Redshift-limited power law 
fits were carried out for a variety of redshifts, of which two particular values z^ax = 2.0, 6.0 
are shown in figures 4 (2B) and 5 (2B-|-PVO) for three values of the evolution index D. The 
case Zfnax = 2 is representative of the redshift inferred by Norris, et al. (1995) based on an 
analysis of BATSE 2B time profiles and brightnesses under the assumption of cosmological 
time-dilation and redshift, without allowance for any possible intrinsic energy-stretching of 
the profiles. The values for z^ax = 6 are characteristic of the maximum redshift inferred 
by Fenimore and Bloom (1995) if there is such an intrinsic energy stretching. We note 
that Mitrofanov, eA al. (1994) have found no evidence for a cosmological time dilation, 
while Norris, et al. (1995) And no strong need for intrinsic energy stretching (see however 
also Fenimore, et al. , 1995). The 2B fits are not bounded from above in the K^Cmax 
plane, due to the lack of very strong bursts in this sample, but especially for the cases 
of strong evolution {D > 3, sources preferentially distant) there is a lower bound to the 
maximum luminosity Umax which is particularly strong for low z^ax^ e.g. for z^ax = 2, 
D = 4: the 3cr limit is Cmax ^ 10^^ s“^ (figure 4c). With the 2B-|-PVO sample, the 
PLZ fits are constrained from above {K^Cmax plane) in all cases considered, as seen in 
figure 5. This is because the strong sources from PVO follow a well defined -3/2 integral 
distribution behavior. For both z^ax = 2 and 6, the 3cr upper limits of K^Cmax are 
~ 10^, 1 — 3 X 10^® s“^ depending on D (figure 5). For % 2 or for z^ax ^ 2, there 
is no restriction against K ~ 1 (i.e. standard candle) or against Cmax ^ Co ^ 10^^ s“^. 
However, the 2B-|-PVO fits constrain K,Cmax from below for % 3 and Zmax ^ 2 (figure 
5b,c). For D = 3 we And K ^ 1, Cmax ^ 10^^ s“^, while for D = A we And K ^ 3, 
Cmax ^ 2 X 10^^ s“^ as the 3(7 lower limits. 

A truncated power law is a highly idealized luminosity function. Nonetheless, it 
is a form commonly found in astrophysics, and its spread can give us some idea about how 
standard is the efflciency of the source in producing y-rays. While it is easy to envisage 
a ’’standard” energy, e.g. ~ Mqc^, it would be more difflcult to envisage converting that 
into y-rays with a well deflned efflciency which is the same in every source. It is thus 
reassuring that, while previous fits have found standard candle models acceptable, and 
limits of a factor ~ 10 can be put on the spread of the distribution containing 90% the 
observed low flux bursts, equally good fits are found for a significant spread {K % 10 and 
up to 300 or 10^ at the Icr or 3a level) in the intrinsic luminosity distribution with a slope 
matching the low flux number distribution. A lower limit on the spread would also be 
interesting, since it might say something about the possible range of masses involved. On 
the other hand, if the luminosity spread was constrained to be large, it might be surprising 
that the turnover of the counts at low P below the —3/2 behavior is not even more gradual 
than what is observed. However, in most cases there is no need for the spread to be very 
large. We And a 3a lower limit K % 3 only for Zmax ^ 2 and % 4, but for Zmax ^ 2 and 
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D ^ 4: there is no significant difference to within 3a between values of 1 ^ K ^ 10^. Such 
conclusions, of course, depend on the accuracy and completeness of the data at low P, and 
one cannot rule out that there may be some as yet undetermined corrections affecting the 
low P counts. A discussion of incompleteness issues and methodological questions is given 
in Loredo and Wasserman, 1995. Additional data may affect any conclusions reached here. 
Here we have confined ourselves to the use of the published 2B data and correction tables. 

The cosmological fits obtained are of good quality (Xred ^ ^ range of 

plausible model assumptions, including both standard candle and truncated power law 
luminosity functions in the brightness or redshift limited cases. The present analysis 
provides a discussion of the statistical significance of such fits including specific density 
evolution and luminosity function parameterizations. The constraints obtained here on 
the evolution parameter and the intrinsic luminosity function spread are stronger than 
in previous analyses due to the inclusion of both BATSE and PVO information, and the 
use of standard deviation measures. However, such constraints do not impose a fine- 
tuning problem, as the allowed parameter space region is of moderate size and includes a 
substantial range of physically reasonable values under the cosmological interpretation. 
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Figure Captions 


Figure 1 : Standard Gandle (SG) cosmological fits for Oq = 1, A = 0 for photon luminosity 
Co (photon/s) and physical density evolution n{z) oc (1 + z)^ (where U = 3 is the 
nonevolving, comoving constant density case). The inner dots are the best fit Xred 
minimum location, with Icr, 2(7, 3(7 contours increasing outwards, a) Top: using the 
BATSE 2B data base; b) Bottom: using the 2B plus PVO information (see text). 

Figure 2 : Power law (PL) luminosity function cosmological fits for Oq = 1, A = 0, 
luminosity index /3 = 15/8 bounded between Cmin and Cmax and physical density 
evolution n{z) oc (1 + z)^ n the brightness limited case, using the 2B sample. The 
abcissa is log^^Q and the ordinate is log^^Q = log^Q^Cmax/Cmin)- The three 

panels a,b,c from top to bottom are for D = 2,3, 4. 

Figure 3 : Power law (PL) luminosity function fits, same parameters as for figure 2, 
brightness limited case, but for the 2B+PVO sample. 

Figure / .• Power law luminosity function (PLZ) fits to the 2B sample for the redshift limited 
case (sources with density evolution and luminosity function as fig. 3, jS = 15/8, but 
only out to finite Zmax)- Left panels: Zmax = 2.0, from top to bottom a),b),c) cases 
D = 2,3,4; right panels: Zmax = 6.0, from top to bottom c),d),e) cases D = 2,3,4. 

Figure 5 : Power law luminosity (PLZ) fits to the 2B+PVO sample, /3 = 15/8, for the 
redshift limited case (other details as for figure 4). Left panels: Zmax = 2.0, from top 
to bottom a),b),c) cases D = 2,3,4; right panels: Zmax = 6.0, from top to bottom 
c),d),e) cases D = 2,3,4. 
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